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Quasiparticle calculations are performed to investigate the electronic band structures of various 
polymorphs of Hf and Zr oxides. The corrections with respect to density-functional-theory results 
are found to depend only weakly on the crystal structure. Based on these bulk calculations as well 
as those for bulk Si, the effect of quasiparticle corrections is also investigated for the band offsets at 
the interface between these oxides and Si assuming that the lineup of the potential at the interface 
is reproduced correctly within density-functional theory. On the one hand, the valence band offsets 
are practically unchanged with a correction of a few tenths of eV. On the other hand, conduction 
band offsets are raised by 1.3-1.5 eV. When applied to existing calculations for the offsets at the 
density-functional-theory level, our quasiparticle corrections provide results in good agreement with 
the experiment. 
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I. INTRODUCTION 

In the microelectronic industry, the continuous quest 
for devices with improved performance and lower power 
consumption has recently stimulated an intense research 
on dielectric materials. Indeed, for over three decades, 
S1O2 has formed the perfect gate dielectric material for 
metal oxide semiconductor field effect transistors (MOS- 
FETs). However, fundamental limits have recently been 
reached that impede further downscaling of MOSFETs 
based on Si02- 1 Candidate materials to substitute the 
latter are transition metal oxides and silicates with an 
high dielectric constant, specifically higher than Si02, 
and are commonly referred to as high-k dielectrics. In 
this framework, Zr02 and HTO2, and more generally Zr 
and Hf compounds, have attracted considerable atten- 
tion,^ hafnium-based microprocessors being now in de- 
velopment or even already in productionP 

Ab-initio calculations can nicely complement the ex- 
perimental work to investigate the properties of these 
novel materials (see, e.g., Refs. 131415ft and to engineer 
the interfaces!^ The method of choice for investigat- 
ing ground-state properties is density-functional theory 
(DFT) that allows to treat quite large systems on the one 
hand, and to obtain reliable results on the other hand. 8 
For the Si/Zr02 and Si/Hf02 in terfaces, various mod- 
els have been explored using DFT) 9 ' 10 ' 11 ' 12 ' 14 ' 15 ' 16 ' It was 
found that due to their analogous electronic structure, 
the two transition metal oxides present a very similar in- 
terfacial bonding. Moreover, there is general agreement 
that the O-terminated interfaces are more stable than 
metal-terminated ones. 

One of the most stringent criteria in the design of 
Si/oxide interfaces is their band offsets (BOs) that con- 
trol the transport properties, and hence the leakage cur- 
rent.^ In particular, both the valence and conduction 
band offsets (VBO and CBO) should be larger than 1 eV 



to obtain a low leakage. DFT relying on local or semilo- 
cal approximations for the exchange-correlation poten- 
tial does not guarantee quantitatively correct BOs since 
the DFT eigenenergies do not corr espond to the quasi- 
particle (QP) energies.^ 8 1 19 ' 20121 1 However, the VBOs 
are often found with an accuracy of a few tenths of eV, 
especially for semiconductor interfaces! 2 ^ Therefore, the 
CBOs can also be predicted using a simple scissor opera- 
tor to correct the band gaps to their experimental values. 
Several works have addressed the BOs at the Si/Zr02 
and Si/ Hf02 interf aces using this scissor-corrected DFT 
schemePSnnEMSI The calculated VBOs for the sta- 
ble insulating O-terminated interfaces of Si/Zr02 and 
Si/Hf02 are around 2.5-3 eV, in reasonable agreement 
with the experiments (2.7-3.4eV)P 24 l 25 l 26 l 27 l 28 l 2t) l 3 ^ 
The scissor-corrected CBOs are about 1.7-2.2 eV, and 
comp are quite well with the experimental values (1.5- 

o c y) j28l29l30l3Tl 

In contrast with DFT, the many-body perturbation 
theory (MBPT) within the GW approximation has 
proven to be a practical and sufficien tly accurate method 
for calculating QP energies j 32 l 33 l 34 l In this method, the 
DFT eigenenergies within the local density approxima- 
tion (LDA) or the generalized gradient approximation 
(GGA) level are corrected perturbatively (QP correc- 
tions) to obtain the QP energies. For semiconductor in- 
terfaces, the QP corrections on the band edges are often 
similar on both sides 3 - and do not substantially affect the 
VBOs, explaining the success of DFTp2l Nevertheless, 
this approach cannot be generalized to other interfaces. 
For example, for the Si/Si02 interface, the difference 
between the DFT and the experimental VBO is larger 
than 1 eV. Recent accurate calculations have shown that 
the QP corrections account for this discrepancy,^ and 
hence they are essential to reproduce quantitatively the 
experimental measurements. For the Si/Zr02 interface, 
a correction of about 1.1 eV has been extracted from GW 
calculations for SP^ and ZrOj^ and used together with 



2 



the experi menta l band gap to correct DFT BOs in sev- 
eral works For the Si/Hf02 interface, the same cor- 
rection as for Si/ZrC>2 has been adopted^ since there 
were no GW calculations available for HTO2. Such an 
assumption seems quite reasonable given the analogous 
electronic structure of ZrC>2 and Hf02- However, for both 
Si/ZrC>2 and Si/HfC>2 interfaces, the VBOs obtained ap- 
plying this correction are too large (and as a consequence 
the C BOs too small) with respect to the available exper- 
imentsPmUS] 

In this work, the origin of this disagreement is dis- 
cussed. QP calculations are performed for the thermo- 
dynamically stable phases (cubic, tetragonal and mono- 
clinic) of HfC>2 and Z1O2 as well for a strained tetrag- 
onal polymorph. In particular, we calculate the QP 
corrections at the top valence and bottom conduction 
bands in order to determine the fundamental band gaps 
and by comparison with Si, the QP correction for the 
Si/oxide band offsets. In contrast with previous QP cal- 
culations, our results show that the VBO remains almost 
unchanged while the CBO is corrected by 1.3-1.5 eV, thus 
explaining the success of the scissor-corrected DFT. The 
paper is organized as follows. In Sec. II, the method- 
ological background is briefly described. Sec. Ill is de- 
voted to the presentation and discussion of our results: 
the DFT geometries and band structures, the QP correc- 
tions to the band gaps, and finally the QP corrections to 
the Si/oxide band offsets. 



II. METHODOLOGICAL BACKGROUND 

The geometries and electronic structures for all the 
systems are computed within the DFT. All the calcula- 
tions are carried out with the ABINIT^3 code within the 
LDA for the exchange-correlation energy functional.^ 
Troullier-Martins^ norm conserving pseudopotentials 
are used which include semicore states in the pseudopo- 
tentials for Zr and Hf (for details see Refs. I3I4[I . The 
wave functions are expanded on a plane-wave basis set 
up to kinetic energy cutoff of 12 Ha for Si and up to 
30 Ha for the HfC>2 and ZrC>2 polymorphs. For all sys- 
tems, the Brillouin zone (BZ) is sampled with a 4x4x4 
Monkhorst-Pack^ grid. 

The QP energies are evaluated using the MBPT within 
the GW approximation. In this approach, the DFT 
eigenenergy E^ FT and wavefunction 4>n FT f° r the n th 
state are used as a zeroth-order approximation for their 
quasiparticle counterparts. Thus, the QP energy E^ F is 
calculated by adding to E^ FT the first-oder perturbation 
correction that comes from replacing the DFT exchange- 
correlation potential v^ FT with the GW self-energy op- 
erator Hgiy: 

e$* = £„ dft + M{ Zn (^ FT \x GW - v» ft \^ ft )}. (1) 

The renormalization factor Z n accounts for the fact that 
EgH?, which is energy dependent, should be evaluated 
at EQ p . The GW self-energy operator Tiqw writes as a 



convolution in frequency space between the one-electron 
Green's function G and the screened Coulomb potential 
W: 

£ GM ,(r,r» = ^ / duj'e iSu ' G(r, r'; u')W(r, r'; uj-u'), 

(2) 

where S is a positive infinitesimal. The explicit expression 
for the Green's function G is: 

G(r ' ' W) " 4- " - ^° FT + « S ^ DFT " ' 

where fi is the chemical potential. The screened Coulomb 
potential is determined as convolution between the in- 
verse of the dielectric function e _1 and the bare Coulomb 
interaction: 



W{r,r';w) = / dr 



r" - r' 



(4) 



The dielectric function e is calculated within the Random 
Phase Approximation. Its dependence on the frequency 
is approximated using the plasmon pole model (PPM) 
proposed by Godby and NeedsP^ For the cubic Zr02, 
we explicitly test the validity of this choice. On the one 
hand, we perform the same calculation with the PPM 
proposed by Hybertsen and Louie.^On the other hand, 
we take into account the full frequency dependence of 
th e diel ectric matrix without resorting to any PPM at 
a llj 4 3 | 44 | m or( j er t discriminate between the two PPMs. 

In our QP calculations, both G [Eq. ((3)] and W 
[Eq. Q] are first evaluated from the DFT eigensolutions 
(which is often referred to as GqWq). Successively, the 
DFT energies are self-consistently replaced in Eq. ^ by 
the corrected values obtained from Eq. (Q (GWq). The 
QP corrections are evaluated only for a few valence and 
conduction bands around the gap, the other energies E^ F 
are extrapolated using a scissor operator. We find that 
2-3 iterations are enough to converge the orbital ener- 
gies up to 0.01-0.02 eV. A systematic study^on several 
bulk systems has recently pointed out GWq as a practical 
and accurate method for evaluating QP energies. Indeed, 
while GqWo provides underestimated energies for almost 
all systems, GWq shows a good agreement with both ex- 
perimental data and full self-consistent GW including a 
Bethe-Salpeter-like vertex correction.^ 

The GW calculations are performed with the yambcJ^ 
code, except for the calculation of cubic Zr02 with the 
PPM of Hybertsen and Louie which is performed with 
the ABINIT code. We carefully study the convergence of 
the QP corrections with the numerical cutoffs: for Si and 
the cubic polymorphs of Hf02 and Zr02, we include 200 
bands in the calculations of the Green's function [Eq. ([3|], 
and 200 bands and 331 reciprocal lattice vectors in the 
calculation of the dielectric matrix^ m Eq. Q . For the 
tetragonal polymorphs (strained and at equilibrium), we 
include 500 bands in the calculations of the Green's func- 
tion, and 300 bands and 735 reciprocal lattice vectors in 
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the calculation of the dielectric matrix. Finally, the mon- 
oclinic polymorphs requires 600 bands for the Green's 
function, and 400 bands and 1177 reciprocal lattice vec- 
tors for the dielectric matrix. With these parameters, we 
estimate an error on the QP energies of about 0.05-0.1 eV 
depending on the system. 



III. RESULTS 
A. DFT geometries and band structures 

For both ZrC>2 and HfC>2, the three thermodynamically 
stable phases (cubic [c], tetragonal [t], and monoclinic 
[m]) are investigated. In addition, a strained [s] form 
of the tetragonal polymorph is also considered in which 
two sides are fixed to a — asi/V% (osi=5-40 A is the 
LDA theoretical lattice constant of Si) while all other 
degrees of freedom are relaxed. This last structure aims 
at simulating the effect of the epitaxial strain on the oxide 
in the MOSFET device. 

The calculated equilibrium parameters describing 
these four geometries are reported in Ta ble Jl] Our re- 
sults aCTeewhhinl-2% with previous LDAEnaHlSl an d 
qqj^MMMMMMM results, as well as with experimen- 
tal dataP^ For the strained polymorphs (s-Zr02 and s- 
HTO2), a contraction is observed along the tetragonal di- 
rection c of about 2% compared to their fully relaxed 
tetragonal analogs (t-Zr02 and t-Hf02). This is a direct 
consequence of fixing the lattice constant a in the basal 
plane (in fact, a is expanded by 7% and 5% in s-Zr02 and 
S-HTO2, respectively). As an indirect consequence, the in- 
ternal parameter d z becomes larger in the strained forms. 
For the latter structure, no direct meaningful compari- 
son is possible with previous works in which epitaxial 
strained tetragonal polymorphs are also considered, ESEU 
since these rely on the GGA value for Si lattice constant. 

The band structure calculated within the LDA are 
reported in brown in Fig. [l] for the cubic, tetragonal, 
strained, and monoclinic forms of Zr02 and Hf02- For 
Zr02, the LDA band structures are in good agreement 
with those presented in Ref. [50] The band structures of 
Hf02 are overall very similar to those of Zr02, with the 
important exception of the cubic phase [Fig. [lja)]. In- 
deed, while c-Zr02 shows an indirect minimum band gap 
from X (top of the valence band) to T (bottom of the con- 
duction band) of 3.4 eV, c-Hf02 has a direct band gap at 
X of 3.5 eV. The indirect X~r gap — that is the minimum 
gap in c-Zr02 — is slightly larger (3.7 eV). 

In the tetragonal phase [Fig. [TJb)], the top valence 
band is almost flat along the TZ and XrM directions. 
It actually presents several maxima (close to X at ^0.20 
in the Xr direction and close to T at ^0.20 in the TM 
direction, in A and Z). The conduction band minimum is 
located at T. The resulting indirect band gap is 4eV for 
both t-Zr02 and t-Hf02, respectively. The effect of the 
strain is to reduce the band gap to 2.5 eV for Zr02, and 
to 3.1 eV for Hf02 [Fig. [ljc)] . The effect is larger for s- 
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TABLE I: Structural parameters of the thermodynamically 
stable phases (cubic [c], tetragonal [t], and monoclinic [m]) 
and of the strained [s] tetragonal structure (see text) of Zr0 2 
and Hf02- The lattice constants (a, b, and c) are in expressed 
in A, while the angle 7 (between a and b) is given in degrees. 
For the tetragonal and strained forms, the internal parameter 
d z is the displacement of the oxygen atoms with respect to 
their ideal cubic position in units of the lattice vector c. For 
the monoclinic polymorph, the internal coordinates for the 
metal (M=Zr or Hf) and the two oxygen (Oi and O2) atoms 
are given in terms of lattice vectors. 



Zr02 since it has the larger mismatch with the Si lattice 
constant. A similar reduction of the band gap resulting 
from t he Si -epitaxial strain was also observed in previous 
works ESEH p or both s-Zr02 and s-Hf02, the valence band 
is almost flat along the XrM direction. The maxima are 
located close to T, while the conduction band minimum 
is at r. 

For the monoclinic phase [Fig. |ljd)], the top of the 
valence band is located at T while the bottom of the 
conduction band is located at B. The indirect minimum 
band gap is 3.7 eV and 3.8 eV for m-Zr02 and m-Hf02, 
respectively. 



B. QP corrections on the band gaps 

For all systems, the QP corrections are reported in Ta- 
ble [II] and the resulting band structures appear in black 
in Fig. [l] The effect of the QP corrections is to lower the 
valence bands and to raise the conduction bands. For 
both Zr02 and Hf02, the QP corrections on the band 
edges depend only weakly on the structure (c, t, s, or 
m). The correction 5E V at the valence band maximum 
(VBM) varies from —0.3 to — 0.5cV in Zr02 and from 
—0.4 to — 0.5eV in Hf02, while the correction 5E C at 
the conduction band minimum (CBM) ranges from 1.4 
to 1.5eV in Zr02 and from 1.5 to 1.7eV in Hf02- More- 
over, the QP corrections are not affecting the location 
of the minima and maxima found in the LDA bandstruc- 
tures. The net effect of the QP corrections is thus to open 
the LDA band gaps by SE g —SE c — 5E V varying from 1.8 
to 2eV for Zr0 2 and from 1.9 to 2.1 eV for Hf0 2 . The 
difference between the Zr02 and Hf02 band gaps found 
in LDA is increased by ^0.2 eV when including the QP 
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FIG. 1: [color online] Theoretical band structure within the LDA (brown circles) and GWo (black circles) along high symmetry 
axis of the BZ for the (a) cubic, (b) tetragonal, (c) strained, and (d) monoclinic phases of ZrC>2 and HfC>2. The GWo band 
structures have been extrapolated from the calculated QP corrections using a linear fit. 
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corrections (from 0.1 to 0.3 eV in the cubic and mono- 
clinic phase, from 0.6 to 0.8 eV in the strained phase), 
except for the tetragonal phase in which the band gaps 
differ by less than 0.1 eV both at the LDA and QP level. 
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TABLE II: For bulk Si, and bulk Zr0 2 and Hf0 2 in the cu- 
bic (c), tetragonal (t), monoclinic (m) phases and strained 
(s) polymorph: the minimum DFT band gap £° FT , the min- 
imum QP band gap E® p at GWo level and the QP corrections 
SE V and 5E C to the VBM and CBM, respectively. 

For ZrC>2, the QP correction to the band gap, which 
is 1.9 eV in the three thermodynamically stable phase, is 
0.4 eV lower than the value previously obtained in Ref.H9l 
for the cubic phase. This variation can be attributed to 
the different PPM (by Hybertsen Louie^l) used by the 
authors in the frequency integration of Eq. Q. Indeed, 
by repeating our calculation for the cubic phase with the 
same PPM, the QP correction to the band gap increases 
up to 2.4 eV in good agreement with Ref. HHl In order 
to discriminate between the two PPMs, the calculation is 
also repeated without resorting to any PPM. The QP cor- 
rection is found to be 2.1 eV, which is 0.2 eV higher than 
the value obtained with the Godby and Needs PPM, and 
0.3 eV lower than the one obtained with the Hybertsen 
and Louie PPM. It can be argued that our QP correc- 
tions to the band gap and hence the resulting QP band 
gaps are probably also underestimated for the other poly- 
morphs of ZrC>2 and for Hf02 due to the use of the PPM. 
Thus, when comparing the results with the experiments 
or other theoretical work, this extra uncertainty of about 
0.2 eV should also be taken into account. 

Our calculated QP band gaps for c-, t-, and m-Zr02, 
also show an overall agreement with those of Ref. 57 re- 
lying on the screened-exchange LDA method. This con- 
firms the validity of this approximation for the calcula- 
tion of band structures. For HfC>2, our result for the 
cubic phase agrees well with the GW band gap reported 
in Ref. 58 , while for the tetragonal and monoclinic phase 
our values are larger by 0.2 eV. 

Experimentally, the optical band gaps determined 
from energy loss and transmiss ion spectros copies range 
from 5.2 to 5.7 eV for Zrq J 31 ' 50 ' 59 ' 60 ^ !!^^ and from 5.3 
to 5.8 eV for Hf0 2 ! 31 | 64 | 65 Mmzl Hence, the agreement be- 
tween our calculations of the fundamental gap and exper- 
imental results is quite reasonable considering the tem- 
perature, excitonic, and impurities effects and the possi- 
ble substrate strain in case of deposited films that should 
be taken into account. For the thermodynamically sta- 
ble phases of ZrC>2 , reflectance vacuum ultraviolet spec- 



troscopy indicate 6.1-7.1 eV, 5.8-6.6 eV, and 5.8-7.1 eV 
for the cubic, tetragonal, and monoclinic phases, respec- 
tively. For the last two, our calculated values fall into 
the range of the experimental estimates. In contrast, 
for the cubic phase, our band gap is much lower even 
when considering the minimum direct band gap (5.5 eV 
at X) and when taking into account the underestima- 
tion by 0.2 eV coming from the PPM. This discrepancy 
may be related to the yttrium used to stabilize the cu- 
bic phase at room temperature and/or the tendency of 
reflectance measurements to overestimate the gap (see 
Ref. \SU\ for a thorough discussion) . For HfC>2 , a "theo- 
retical" band gap of 6.7 eV^ was proposed for a film de- 
posited on a SiO^Ny/p-Si substrate (monoclinic phase) 
by comparing direct/inverse photoemission spectroscopy 
with DFT density-of-states calculations, arguing that the 
reduction of 0.9 eV with respect to the experimental band 
gap (5.9 eV) should be attributed to defects tail states. 
Our value compares better with the band gap determined 
directly from the experiment. 

C. QP corrections on the Si/oxide band offsets 

In the DFT approach, the VBO and CBO are conve- 
niently split into two terms: 

VBO = A£ DFT + AV (5) 
CBO = AE? FT + AV (6) 

The first term A£'° FT [resp. AE® FT ] on the right-hand 
side of Eq. ^ [resp. Eq. ([6])] is referred to as the band- 
structure contribution. It is defined as the difference be- 
tween the VBM [resp. the CBM] relative to the average 
of the electrostatic potential in each material. These are 
obtained from two independent standard bulk calcula- 
tions on the two interface materials. The second term 
AV^, called the lineup of the average of the electrostatic 
potential across the interface, accounts for all the intrin- 
sic interface effects. It is determined from a supercell 
calculation with a model interface. 

Despite the DFT limitations in finding accurate 
eigenenergies, the VBOs are often obtained with a very 
good precision, in particular for semiconductors.^ This 
has opened an indirect route to compute the CBOs 
through the experimental band gaps using: 

CBO = AE c g xp + VBO (7) 

where AEg Xp is the difference between the experimental 
values of the band gap of each material. Note that this 
equation is equivalent to applying a scissor correction to 
the conduction bands on both sides of the interface, as 
can be seen by inserting Eq. ^ in Eq. 

CBO = A£ C DFT + AV+ (AE c g xp - A£° FT ) (8) 

and comparing with Eq. ([6]). 

As discussed in the introduction, this scissor-corrected 
DFT scheme has been used in several studies of the 
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Si/Zr0 2 and Si/HfOaP 121111111 For the stable insulat- 
ing O-terminated interfaces of Si/ZrC>2 and Si/HfD2, 
the VBOs calculated are found to range from 2.5 
to 3 eV, in good agr eement with the e xperimen- 
tal findings (2.7-3.4 eV).^ 24 ' 25 ' 26 ' 27 ' 28 ' 29 ' ™ Adopting 

&E° xp =4.7eV (^r p = !- leV for Si and 5 - 8eV for Zr °2 
and Hf02) in Eq. Q, the scissor-corrected CBOs lie thus 
between 1.7 and 2.2 eV. This is also in good agreement 
with experiments (1.5-2.0 eV) J 28 ' 29 ' 30 ' 31 ' 

In the QP framework, it was often assumed^ and it has 
recently been proverPS that the lineup of the potential 
AV is already well-described within DFT. So that, only 
the band-structure contribution is modified: 

VBO = AE$ P + AV = AE^ FT + A(SE V ) + AV (9) 
CBO = A£ C QP + AV = A£ C DFT + A(5E C ) + AV (10) 

where 6E V = E$ p - E» FT (resp. 5E C = E? p - £ C DFT ) 
is the quasiparticle correction at the VBM (resp. CBM) 
and A(SE V ) [resp. A(8E C )] is the corresponding differ- 
ence between the two materials. 

For the Si/Zr0 2 and Si/HfC>2 interfaces, no specific 
GW study exists as such. However, Fiorcntini et al™ 
evaluated the QP effects on the VBO of Si/ZrC>2 combin- 
ing the value for tetragonal Zr0 2 (SE V =— 1.23 eV) from 
Ref.[49]with the value for Si from Ref.|23(6£;. u =-0.15 eV. 
This results in a total correction of A(SE v )=l.08eV on 
the VBO in Eq. This value has been used in several 
other worksp2H3 even for Si/Hf02 interfaces. For both 
Si/Zr02 and Si/Hf02 interfaces, the VBOs obtained in 
this way were found to be too large (and as a consequence 
the C BOs too small) with respect to the experimental 

vaiueBpannni 

For the oxides, our QP corrections to the DFT va- 
lence band SE V vary from —0.3 to — 0.5eV; while for 
the conductions bands, SE C ranges from 1.4 to 1.7 eV. 
For Si, our QP corrections, which are reported in Ta- 
ble [II] lead to a band gap that agrees well with previous 
theoretical works (e.g. see Ref. |4"5| . and with the ex- 
perimental value™ The total QP correction to the gap 
5E g of about 0.7 eV, comes mostly from the downshift of 
valence band state (5E V =— 0.6eV). This 6E V value is al- 
most the same as that found for the polymorphs of both 
Zr02 and Hf02- Therefore, the QP correction on the 
VBOs are only of 0. 1-0.2 eV and the correction on the 
CBOs is about 1.3-1.5 eV. This explains why previous 
studies based on scissor-corrected DFT were in such a 
good agreement with experi mental r esults. 

Turning to previous works^EsEll that accounted for 
QP corrections to the BOs of Si/Zr02 and Si/Hf02 inter- 
faces using values from prior GW calculations! 3 ^* 4 ^ their 
disagreement with experiments can be explained as fol- 
lows. On the one hand, the QP corrections in the oxide 
and Si are not consistent since a different approximation 
has been used in Eq. Q for the dielectric function e. 
The calculations for Si use a model dielectric function,^ 
while the calculations for Zr02 use the Random Phase 
Approximation. In particular, the value obtained for Si 
5E V =— 0.2eV is lower in absolute value with respect to 



the value found from the Random Phase Approximation 
(5E V =— 0.6eV), and hence this inconsistency artificially 
increases the QP correction on the VBO. On the other 
hand, our QP results (in particular, those for c- and t- 
Zr02) differ from those obtained in Ref. l4"9l where it was 
found that the total QP correction SE g =2.3eV resulted 
from lowering the valence bands by about 1.3 eV, and 
raising the conduction band by about 1.1 eV. The differ- 
ence with our results is due again to the different PPM 
used. Indeed, repeating the calculations with the Hy- 
bertsen and Louie PPM for the cubic phase, we found 
a correction 6E V of — 1.1 eV for the valence and SE C of 
+1.3 eV for the conduction in agreement with Ref. 
When performing the same calculation without resorting 
to any PPM, we obtain &E„=-0.7eV and 5E c =lAeV, 
meaning that SE V is 0.2 eV too low for the Godby and 
Needs PPM, and 0.5 eV too high for the Hybertsen and 
Louie PPM. For the conduction bands, 5E C differs by 
less than 0.1 eV with the former, while it is 0.3 eV too 
low with the latter. 

As final remark, we stress that rigorously, the QP 
corrections on the band offsets should be calculated us- 
ing the same pseudopotential and the same exchange- 
correlation approximation as for the interface calcula- 
tions. Indeed, while the GW band gap has been demon- 
strated to be quite insensitive to the starting point, 
this is not true for the QP corrections that reflect — 
as a consequence — the differences in the choice of the 
pseudopotential and the exchange-correlation approxi- 
mation. Therefore, we would recommend to calculate — 
when possible — QP corrections on DFT band offsets us- 
ing the same pseudopotentials, and exchange-correlation 
approximation as for the interface calculation and using 
the same PPM for both materials. 



IV. SUMMARY 

The electronic properties of Zr02 and Hf02 poly- 
morphs and their interface with Si have been investigated 
using GW calculations. The QP corrections are found 
to be very similar for the two oxides consistently with 
their analogous band structure, and depend only slightly 
on the crystalline structure. While, at the DFT level, 
the epitaxial strain was found to dramatically shrink the 
band gap (especially for Zr02 for which the lattice pa- 
rameter mismatch with Si is larger), the QP corrections 
depend only slightly on the strain. When considering 
the interface Si/oxide, the QP corrections to the VBOs 
were calculated to be very small (a few tenths of eV) by 
cancellation of the corrections on the valence band max- 
imum of the Si and those of oxides. On the other hand, 
the correction was found to be of the order of 1.5 eV for 
the CBOs. These results disagree with the correc tion on 
the VBOs of more than 1 eV used in the literature] 14 * 15 * 16 ! 
which was extracted from existing GW calculations! 35 * 49 ^ 
We have traced back the differences with our results to 
the difference in the PPM used for the Zr02, and to 
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the inconsistence in the level of approximation for the 
screened interaction W. Our results combined with the 
DFT band offsets available from the literature for differ- 
ent interfacial bonding models provide values in a good 
agreement with the experiment. 
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